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Abstract 

In this paper, with particular emphasis on the N-body problem, we focus our attention on the possibility 
of a synthetic characterization of the qualitative properties of generic dynamical systems. The goal of this line 
of research is very ambitious, though up to now many of the attempts have been discouraging. Nevertheless 
we shouldn't forget the theoretical as well practical relevance held by a possible successful attempt. Indeed, if 
only it would be concevaible to single out a synthetic indicator of (in)stability, we will be able to avoid all the 
consuming computations needed to empirically discover the nature of a particular orbit, perhaps sensibly different 
with respect to another one very near to it. As it is known, this search is tightly related to another basic question: 
the foundations of Classical Statistical Mechanics. So, we reconsider some basic issues concerning the (argued) 
relationships existing between the dynamical behaviour of many degrees of freedom (mdf) Hamiltonian systems 
and the possibility of a statistical description of their macroscopic features. The analysis is carried out in the 
framework of the geometrical description of dynamics (shortly, Geometrodynamical approach, GDA), which has 



been shown, |Cipriani 1993, Pettini 1993, Pi Bari 1996, Cipriani et al. 1996, Pi Bari et al. 1997 1, to be able to 



shed lights on some otherwise hidden connections between the qualitative behaviour of dynamical systems and the 
geometric properties of the underlying manifold. We show how the geometric description of the dynamics allows 
to get meaningful insights on the features of the manifold where the motion take place and on the relationships 
existing among these and the stability properties of the dynamics, which in turn help to interpret the evolutionary 
processes leading to some kind of metastable quasi-equilibrium states. 

Nonetheless, we point out that a careful investigation is needed in order to single out the sources of instability on 
the dynamics and also the conditions justifying a statistical description. Some of the claims existing in the literature 
have here been reinterpreted, and few of them also corrected to a little extent. The conclusions of previous studies 
receive here a substantial confirmation, reinforced by means of both a deeper analytical examination and numerical 
simulations suited to the goal. These latter confirm well known outcomes of the pioneering investigations and 
result in full agreement with the analytical predictions for what concerns the instability of motion in generic 
classical Hamiltonian systems. They give further support to the belief about the existence of a hierarchy of time- 
scales accompanying the evolution of the system towards a succession of ever more detailed local (quasi) equilibrium 
states. These studies confirm the absence of a direct (i.e. trivial) relationship between average curvature of the 
manifold (or even frequency of occurrence of negative values) and dynamical instability. Moreover they clearly 
single out the full responsibility of fluctuations in geometrical quantities in driving a system towards Chaos, whose 
onset do not show any correlation with the average values of the curvatures, being these almost always positive. 
These results highlight a strong dependence on the interaction potential of the connections between qualitative 
dynamical behaviour and geometric features. These latter lead also to intriguing hints on the statistical properties 
of generic many dimensional hamiltonian systems. It is also presented a simple analytical derivation of the scaling 
law behaviours, with respect to the global parameters of the system, as the specific energy or the number of degrees of 
freedom, of the relevant geometrodynamical quantities entering the determination of the stability properties. From 
these estimates, it emerges once more the strong peculiarity of Newtonian gravitational interaction with respect 
to all the topics here addressed, namely, the absence of any threshold, in energy or number of particles (insofar 
N2> 1) or whatsoever parameter, distinguishing among different regimes of Chaos; the peculiar character of most 
geometrical quantities entering in the determination of dynamical instability, which is clearly related to the previous 
point and that it is also enlightening on the physical sources of the statistical features of N-body self gravitating 
systems, contrasted with those whose potential is stable and tempered, in particular with respect to the issue of 
ergodicity time. We will also address another elementary issue, related to a careful analysis of real gravitational 
N-body systems, which nevertheless provides the opportunity to gain relevant conceptual improvements, able to 
cope with the peculiarities mentioned above. 
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Introduction 



The discipline of Deterministic Chaos, in the last decades, has known an expanding phase hardly comparable with 
other fields of mathematical physics. The increased evidence of the ubiquity of instability has spurred an improved 
understanding of nonlinear dynamics. However, at present, as usual for quickly expanding areas of research, after 
the first phase of spontaneous growth, the demand for a more systematic settlement, together with a knowledge of 
the chaotic phenomena going beyond the one of semi-phcnomcnological nature, is also rising. The signature of Chaos 
in realistic models of physical systems is indeed often phenomenological in character, as the occurrence of instability 
is usually recognized a posteriori, looking at its consequences on unpredictability through the study of the evolution 
of a generic perturbation (which represents, from a physical viewpoint, the unavoidable uncertainty on the initial 
conditions) to a given reference trajectory. However, recently there has been a renewed interest in looking for, at 
least qualitative, synthetic indicators of dynamical instability in realistic models of physical systems, able to provide 
an a priori criterion. As far as we know, the first explicit concrete attempt in this direction dates back to Toda's 
paper, Toda 1974 1, whose goal was to find a concise characterization of the sources of global instability related to 
local criteria. As it has been shown (see, e.g., [Benettin et al. 1977, Casati 1975 ), the relationship argued by Toda 
was somewhat too naive, and the link between local and global instabilities turns out to be by far more complex than 
sought. 

Nevertheless, mainly in the field of celestial mechanics, where the dynamical systems under study possess few degrees 
of freedom {few, here, means of order ten, or less), a number of attempts to single out synthetic indicators of instability 
(although qualitative), have been put forward, looking at the intermingled structure of deformed or destroyed invariant 
surfaces in phase space. This effort has led to a deeper understanding of the conditions responsible for the onset of 
Chaos in dynamical systems of low dimensionality and, in a sense, not too far from integrable ones. 
At the other end, very far from integrability, within the sphere of Ergodic Theory, with reference to Statistical Me- 
chanics (although not so close to the latter as the very origin and motivations of the former should let to expect), 
and consequently mainly for mdf systems, some results have been obtained in the case of abstract dynamical systems, 
[ Anosov 1967 , ginai 1991 1, but they seem to be of little utility to gain some insights on the behaviour of realistic 
physical models. Whereas, being mainly focused on conceptual issues on the basic postulates of Analytical and Sta- 
tistical Mechanics, and for this reason relevant to the present work, the last decades knew a deeper reconsideration of 
the issue about the origin of unpredictability and of the problem related to the approach to equilibrium for generic 
mdf Hamiltonian systems (see, e.g., [ Galgani 1988 | and references therein). This problem arose after the first nu- 
merical experiments performed on a computer, when Fermi, Pasta and Ulam, |Fermi et al. 196f:], tried to follow 
numerically the dynamics of a chain of weakly coupled nonlinear oscillators, in the hope of observing the theoretically 
expected process of equipartition of energy among normal modes. Those experiments, because of their surprising 
results, stimulated a lively debate, and supported the believe that the semi-heuristic argumentations of Boltzmann 
himself and Jeans about the exponentially long equipartition times amongst high and low frequency modes possess a 
very deep physical meaning^. 

Among the attempts towards intrinsic (or synthetic) criteria to single out the presece of Chaos in models of concrete 
physical systems, besides the Toda's trial already mentioned, we recall the Krylov's Ph.D. thesis, [Krylov 1979|, whose 
intuitions had important influence on the whole line of research this work belongs to. 

Nevertheless, the founding arguments of the investigation of the stability of motions by means of a geometrization 
of the dynamics date back to the paper by Synge, |Synge 1926 1, which also partially summarizes previous works of 
Ricci and Levi-Civita. Although completely unaware of the relevance of exponential instability (of course, in 1926!), 
he there gave a very clear exposition of the power of the Geometrodynamics in gaining intrinsic information on the 
qualitative behaviour of generic holonomic dynamical systems. Other authors, jEisenhart 192£], since then, added 
their own contribution to the development of the subject, extending the range of applicability of the method. In the 
recent years a new deal of interest has been raised, stimulated by the Krylov's work cited above, (unfortunately diffused 
worldwide more than two decades after its appearance) and from the mess of papers produced by the abstract ergodic 
theory (see, e.g., [ 3inai 1991 1). But only very recently the approach has received a self-consistent settlement able to 
cope with the occurrence of instability in mdf hamiltonian systems (see [ Pettini 1993 , Cipriani 1993 |) and to give a 
key to understand the links between local properties of the manifold and the asymptotic qualitative features of the 
dynamics. In a recent paper, [Di Bari et al. 1997 1, we extended further the Geometrodynamical approach, to handle 
a wider class of dynamical systems, even with peculiar lagrangian structure, embedding the dynamics in an ensemble 
of manifolds more general than riemannian: the Finsler ones. This framework enable us to cope with dynamical 

^For a very interesting account on these pioneering and, for a long time, rather disregarded intuitions, see again [Glalgani 1988 1. 
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systems otherwise not suitable for a geometrization within the riemannian geometry, [Di Bari and Cipriani 1997a |, 



to get more insights on the nature of instability, to locate the relationships with the geometrical structure of the 



manifold, | Cipriani and Di Bari 1997c , and to set up a formalism with a built-in explicit invariance with respect to 



any reparametrization of the time variable (see also |Di Bari and Cipriani 1997b ). The GDA has proved to give many 



interesting insights on the analysis of qualitative properties of dynamical systems (DS) of interest in various fields 
of Statistical Mechanics, as well in Celestial Mechanics and Cosmology; we will address here a problem that, also 
looked from a number criterion, seats just in the middle between Analytical (log ^ 1) and Statistical (log TV ~ 23) 
Mechanics. The main concern of the present paper is indeed towards a clarification of the links between geometry and 
dynamics on one side, and their implications on the statistical behaviour of N-body systems, when 6 <, log <, 13, 
on the other. In the sequel, we will focus on the properties of the gravitational N-body system and, using the Fermi- 
Pasta-Ulam (FPU) chain as a normal reference system, belonging to the class of mdf systems for which a Statistical 
Mechanical description is well founded, compares their dynamical, geometrical and statistical features, enlightening 
from a perhaps nouvel perspective the physical and mathematical peculiarities of Gravity, showing also how to treat 
successfully some of them. 

Instability and Statistical description. 

The relevance of dynamical instability and ergodicitj^ to Statistical Mechanics (SM) has been sometimes questioned, 



[ Farquhar 1964 1 , but it is by now generally accepted that the unpredictability implied by Chaos is a natural ingredient 
to justify the statistical description of mdf dynamical systems. As stressed in the introduction, however, there is no 
rigorous, or even convincing, prove of any direct relationship between the instability times, as can be derived from 
dynamics (Lyapunov exponents, Kolomogorov-Sinai entropy,...) and the time-scales related to a Statistical Mechanical 
treatment, i.e., those belonging under the name of relaxation times. Nevertheless some qualitative relationships exist 
among them for some interesting and historical mdf systems, as the FPU chain, for which, although differing each other 
by orders of magnitude, the Lyapunov and the SM relaxation times display the same transition between two different 



regimes of weak and strong chaoticity, |Pettini and Cerruti-Sola 1991]. Recently (see [Morbidelli and Froeschle 1996 



and references therein), also for dynamical systems with few degrees of freedom and of particular interest in Celestial 
Mechanics, some debate raised about such a relation among Lyapunov exponents and the time scale of diffusion in 
phase space. 

A thorough treatment of the multiple connection between dynamical instability and Statistical Mechanics is out 



of place here, [Cipriani 1993, Cipriani and Di Bari 1997b . Instead, we would like to stress that the suggestions of 



Boltzmann and Jeans, reminded in the previous section, received, in the framework of Analytical Dynamics, a strong 
support from the Nekoroshev theorem on the exponentially long time scales of quasi- conservation of actions for 



generic non integrable systems, jBenettin 198S|. This is relevant to the present discussion in giving an explanation 
to the existence of different regimes of stochasticity in generic Hamiltonian systems of interest in Solid State Physics 
as well in Astrophysics, and also in demonstrating once more the peculiarity of Newtonian Gravity. Indeed, for the 
N-body system governed by gravitational interaction, we will support the claim according to which simply do not exist 
a regime of quasi integr ability, as this system will be chaotic for almost all generic initial conditions^. The absence 
of any integrable limit, makes it hard to find out any threshold, e*, for the magnitude of the perturbation, e, below 
which the actions are quasi conserved for time intervals exponentially long in some power of (e^/e); simply, we don't 
know the state for which e = 0! 

Behind the long standing problem of the justification of the basic assumptions at the ground of the Statistical 
description, fMa 1988| , for generic dynamical systems, lies a slowly evolving set of ideas, all referred to as the Ergodic 
hypotesi^ The GDA, applied to N-body gravitational system, discloses evidence of the differences amongst some of 
the interpretations given of it in order to justify the approach to equilibrium, and, at the same time, reinforces the 
claim that the processes leading a DS to the hierarchy of intermediate partial equilibria, proceed on time scales whose 

^We do not intend here neither to confuse nor to assimilate the two properties. We would only point out how even some very basic 
questions concerning the links between Analytical and Statistical Mechanics are still awaiting for an universally accepted, although not 
rigorously proved, answer. 

•^We are supposing obviously that N 1, and we speak about chaos in a not too strict sense, as the long range behaviour of the 
interaction, more than its short range divergencies, causes the phase and configuration spaces of the system to be non compact. 

■^The ensembJe, often very inhomogeneus, of concepts indicated under this name has sometimes been presented in a somewhat misleading 
fashion, leading, from time to time, either to evidently unreliable assumptions or to trivial assertions. This has been an undeserved 
treatment to Boltzmann efforts to find a strict formulation for a very deep and notwithstanding ineffable physical intuition. A thorough 
and enlightening conceptual re-examination of the issues related to the Boltzmann ideas on the importance of a physically reasonable 
formulation of the Ergodic hypotesis to justify the use of Statistical a pprnar.h in Mech anics has been carried out recently in a scries of 
paper by G.Gallavotti; for an exhaustive account, see the monography JGallavotti 1995 
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spectrum is written in the interaction potential governing its dynamics, and can differ by orders of magnitude passing 
from a DS to another. We refer to [Cipriani 1993, Cipriani and Di Bari 1997b| for a less synthetic discussion of 
these issues, and concentrate here mostly on the dynamics of N-body systems and on the insights on their qualitative 
properties that can be obtained in the framework of GDA. 



Geometrodynamics of evolutionary processes. 

The problem of the approach to equilibrium in physical systems with a large number of degrees of freedom is, essentially, 
a problem of time-scales. The fundamental question is: what time-scale is required such that the system finds itself 
in a state in which some or all memory of the initial one is lost? A theory of non-equilibrium processes that should 
be able to give quantitative answers to this question would also provide the hints to solve most of the open issues 
referring to the connection between Analytical and Statistical Mechanics. This problem has been investigated since 
the very birth of SM with different methods; here we will exploit the Geometrical approach. 

In the following sections we recall briefly the method and then introduce the main Geometrodynamical Indicators 
(GDI's), whose behaviour determine the qualitative (global) evolution of large N-body systems. 

In order to highlight the singular features of Gravity, we apply the method to two potentials, representative of the 
classes of short and long range interactions respectively. The evaluation of the characteristic time-scales for them is 
firstly derived on the basis of analytical estimates of the quantities entering the GDI's. On these grounds some of the 
existing results are reviewed and reinterpreted, giving an (hopefully) coherent scenario. After that, the reliability of 
the claims based on a trivial extension of methods and concepts borrowed from Ergodic Theory are critically examined 
and discussed. Still within a semi-analytical approach, we show why that extension results, at least, too much nai've. 
Using numerical simulations suitably tailored to the goal, we then show that the analytical estimates correctly describe 
the behaviour of GDI's, and so that the mechanisms driving to Chaos in mdf Hamiltonian systems differ completely 
from those occurring within the Ergodic theory of abstract DS's. 

We conclude comparing the behaviour of the GDI's for the mathematical Newtonian potential, on a side, with that of 
the corresponding quantities related to the physical gravitational N-body problem and to the FPU chain, on the other. 
This allows us to shed some lights on the way a stellar systems could attain its apparent equilibrium state, in spite of 
the celebrated theorems of rigorous Statistical Mechanics on non-existence of a maximum entropy state for unscreened 
Coulombic interactions, which do not possess the property of stability, [ Ruelle 1988| . The search for the conditions 



needed in order to guarantee the equivalence between time and phase averages of geometric quantities, leads also to 
the conclusion that the peculiarities of Newtonian gravity, reflects themselves in a singularity on the literal ergodicity 
time, singularity which disappears when a realistic interaction, that nevertheless do not modify the very nature of 
Gravity, is considered, irrespective to the details of the care. 

Though there are several geometrizations of Dynamics, in the following, we exploit mostly that based on the so-called 
Jacobi metric, which results from a straightforward application of the least action principle in the form given by 



Maupertuis, see, e.g., Goldstein 1980, Arnold 1980 1. 



The Jacobi Geometrodynamics. 

In the framework of a definitely non-perturbative treatment of the Hamiltonian Chaos, the differential geometrical 
approach stands in the pathway that, started by Krylov with the aim to comprehend the relaxation processes in realistic 
physical systems, has grown very long in the applications to abstract dynamical systems, acquiring mathematical 
strictness but only recently broadening the physical interest. The Geometrodynamical approach to a Hamiltonian 
system reduces it to a geodesic flow over a manifold, M, on which is defined a suitable metric. In the case of the Jacobi 
geometrization, the manifold is a riemannian one, equipped with the (conformally flat) metric, g = {5ab}, whose line 
element is 

dsf =^ gabdq'^dq'' = [E - Uiq)]7^abdq''dq' = 2W^ dt\ a, 5 = 1, . . . , ; (1) 

where E is the total conserved energy, U{q) is the potential energy, depending on the coordinates q = on the 
configuration manifold of the system (with Af degrees of freedom), 'W{q) [E — Z^(q)] is the conformal factor, 
whose magnitude numerically equals given the kinetic energy T, and t is the newtonian time. The metric of the 
AA-dimensional physical space of the system is thus rjab, i-c., 2T = rjabq°'q^ , which reduces to rjab = ^ab when the space 
is euclidean and cartesian coordinates are employed; as usual, we denote with a dot differentiation with respect to 
time. Since we are considering N-body Hamiltonian systems, we have A/" = N • c?, where d is the dimensionality of the 
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system. So, for the self-gravitating system we have d — 3, whereas = 1 for the FPU chain. 

In the language of GDA, the Af second order Lagrangian equations of motion (or the 2Af first order Hamiltonian ones) 
are replaced by the Af geodesies equations 

Z_,_ + rj.„V = „, = i ^-i p) 

where u'^ ~ dq"" /ds, V /ds is the total (covariant) derivative along the flow, the T^^ are the Christoffel symbols of the 
metric g and the summation convention is understood. 

For the metric defined in (Q), namely gab = Vy(q)77a6, and using cartesian coordinates, the geodesic equations read^ 



d'^q" 1 
~ds^ ^ 2W 



'dWdq'^dq'' ^,BW 

2 g 

dq'^ ds ds dq'^' 



0, , (a=l,...,AA) (3) 



and give the trajectories on the configuration manifold in terms of the affine parameter s (conciding with the Maupertuis 
action). The iiewtonian time law of percurrence of the trajectory is obtained via the relation ds = \^yV{q) dt, 
exploiting which one obtain the familiar equations of motion, g° + U"' = 0. 

Once rephrased the dynamics as a geodesic flow, we are left with the determination of the question of stability of 
geodesic paths on the Jacobi manifold. Within the framework of Hamiltonian dynamics this issue is addressed using 
the tangent dynamics equations, which determine the evolution of a small deviation vector, D — (^, ^), in phase-space, 
describing a perturbation to a given trajectory. 

Although the equations for the trajectories coincide with geodesies ones once rephrased in terms of the same parameter, 
the equations for the displacements, which are those relevant to the issue of stability of motion, differ in general. We 



refer to |Di Bari and Cipriani 1997b for a discussion of the problem of equivalence among different linearized equations 
for displacements and here simply remind that all geodesies have, by definition, unit velocity vector, so, in a sense, 
the deviation between geodesies gives the distance between paths and not between points on them. Another point to 
be remarked concerns the apparent self consistency of the geometrical description, which allow for a built-in distance, 
lacking in the usual definition for the norm of a vector in phase space, where an Euclidean structure is imposed by 
hand. 



Geometric description of instability. 



As it is well known, in the framework of Riemannian geometry^ the evolution of a perturbation, Sq, to a geodesic is 
described by the Jacobi-Levi-Civita (JLC) equation for geodesic spread: 

V fVSq" 



ds V ds 



n\Sq'^o, 



(4) 



where it has been defined the stability tensor ?i°ci defined along any geodesic and related to the Riemann curvature 
tensor, R"'bcd, associated to gab by 



ri c — n bcdU U , 

where u is the unit tangent vector to the geodesic. In the case of Jacobi metric 7i reads 

(gradW)2 



(5) 



4>V2 [~dij ^ 2W ds^ 



4>V3 



-U Uc + 



1 

2W 
3 

'4W2 



(6) 



fdW\ 



which clearly depends also on the position, q e M. Here, as usual, VV,a = dW/dq"", the summation convention is 
adopted, and the Euclidean A/'-dimensional gradient operator acting on a function / defined on M has been indicated 
explicitly as [grad/], rather than V/, to avoid confusion with the already defined total derivative along the flow V /ds. 

^Since we restrict to consider mostly Jacobi geomctrization, we will from now on omit the subscript in dsj whenever there is no risk of 
ambiguity, writing simply ds. 



"And also in more general differential manifolds, see [Di Bari 1996, Di Bari et al. 1997 1 
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Equation (Q) contains all the information about the evolution of a congruence of geodesies emanating within 
an initial distance z = ||i5q|| = (5ab 5q°'5q^Y^'^ from the reference one. The problem for a system with mdf is 
that the Riemann tensor contains ©(A/""*) components, and the evaluation of the stability tensor H, is an huge 
task. Nevertheless, just when the dimensionality of the ambient space is high, some assumptions can be made 
on its global properties, averaging in a suitable way. This procedure has been adopted in the past (see, e.g.. 



[ Pettini 1993 , Cipriani 1993 and references therein) with success in describing the behaviour of mdf systems of interest 
in Statistical Mechanics. Indeed, if we are interested in the stability properties of the flow, the relevant quantity is 
the magnitude of the perturbation, measured within the GDA by its norm. Indeed, if = zv, where v is & vector of 
the unitary tangent space at q, TqM, the norm z evolves according to: 



ds2 



(7) 



which is still an exact equation derived directly from eq.(^), [Cipriani 199S |. However, eq.(|^) contains still the full 
stability tensor, and the task has so been made easier only partially, reducing to a single one the number of equations 
to be integrated. To close this equation some assumption is needed. 



Geometric Indicators of Instability. 

We refer to [ Cipriani 1993 Pettini 1993| and subsequent works for the description of the details of the averaging 
procedure, and we want here recall only some critical points which deserve more cautious treatment than the one 
given sometimes, |Gurzadyan and Savvidy 1986 , Kandrup 199C|. Indeed, one intuitive approximation that help to 
partially get rid of the huge amount of computational task in the evaluation of the stability tensor Ti. is suggested 
by the physical reasoning that we are interested in the evolution of a generic small perturbation, however oriented 
with respect to the given geodesic. Actually, the 7t tensor contains all the informations about the local evolution of 
any perturbation to the system, local in the sense that it describes the behaviour of the vectors of tangent space in 
a neighbourhood of the state (q, u). So, a first averaged equation can be derived from eq.(|^) letting Sq to have a 



random orientation with respect to the given Bow. So, as we have by definition that ||f | 



9abV V 



1, being the 



Jacobi metric conformally flat, this equation would be consistent with a randomly chosen orientation for u if we put 
(y°-ii^) = g°-''/J\f. Nevertheless, while practically irrelevant when iV ^ 1, we should make here a comment about the 
correct way of averaging the deviation over directions. Obviously, given a geodesic, there are (A/" — 1) independent 
orientations orthogonal to u, and we can imagine also a perturbation with a non vanishing component along the 
geodesic itself. This straight argumentation leads to the Af~^ factor in the average above. Yet, it follows directly from 
the very definition of the stability tensor, H, i.e., from the simmetry properties of Riemann tensor, that at least one 
amongst its local eigenvalues vanishes, being that associated with the direction along the flow: 



rCc = 0, 



V l = Pu. 



(8) 



with P{s) any scalar quantity. So, a deviation along the geodesic cannot evolve more than linearly in the s-parameter. 



and consequently do not contribute to exponential instability. Moreover, it is known, e.g., |Katok and Hassclblatt 1995 
§17.6], that the evolution of any parallel component decouples from that of the remaining normal one Zj_; then, as 
interested in the possibly exponential growth of the deviation, we are left with only J\f ~ 1 independent orientations 
of v, once fixed the geodesic. Said in other words, we consider only perturbations orthogonal to the geodesic, so that 



gabiy^u' = 



which leads to 



ah 



Within this reasonable and self consistent approximation, equation (^ reduces to the averaged one: 
d^z , , , 

-j^ + Ku(q) 2 = 

where 

A/"-! 



fcu = fcu(q)= AA^^^'^^AA^Tr-H = -^^if(2)(u,e^); 

A=l 



(9) 



(10) 



(11) 
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from now on, as above, the dependence on q G M of all geometric quantities will be understood whenever no confusion 
can arise. 



We now sketch the line along which it can be shown, | Cipriani 1993 , that this quantity equals the (normalized) Ricci 
curvature in the direction of the flow: 

, . dof Ric(u) 

- JJ^-^Y ' (12) 
being the Ricci curvature along any direction indicated by the unit vector ba G TqM, defined as 

Ric(e^) 1^'' Rbce'Ae'x , (13) 
where Rac R''abc = g^'^Rabcd IS the Ricci tensor. 

The quantities i^(^)(u, e^) in eq.(pl]) are the sectional curvatures in the (TV — f) independent planes spanned by u 
and the {e^}, {A — I, . . . ,N ~ I), (e^, e^) = 5ab, (©a, u) = 0. In general, given a riemannian manifold, (M, g), the 
definition of sectional curvature in the plane spanned by two non parallel vectors, (x,y) € TqM is defined as 

l|x|r lly|l - l|x-y|r 

where the norm and the scalar product are those induced by g. 

Without loss of generality, [ 3ynge and Schild 1978| , we can suppose x _L y; so, given an A/'-dimensional manifold, by 
simmetry it follows that there are Af{Af — l)/2 independent sectional curvatures. Further, given an orthonormal basis 
in M, {e^}, {A — 0, ... ,Af — 1), we indicate these sectional curvatures as 

Ki^) if(2)(e^,ej ^ i?.,,, e>^^e>^, ^ 4^), (B ^ A). (15) 

We preliminarly show that the Ricci curvature along the direction defined above, eq.(|l3|) is also equal to 

Trlt'^Ti^.Y^'eiei (16) 

B=0 

where we defined the M tensors jt \d Rabcd s-iid the prime stands to indicate that the sum is over B ^ A. 

Indeed, in order to see that the r.h.s. of eq.(^ actually equals the trace of the corresponding tensor, we observe that, 
by definition, 

AA\ \ - ' b d T3 a c \ " ' b d 

and, using the anti-simmetry properties of the Riemann tensor, we can include in the summation also the term B — A^ 
which vanishes identically. 

Being {e^} an orthonormal basis, we then have 

Y.«-9'", (17) 
SO that: 

Rabcd e>^^ = '^''bd g'" = li^'^b - TrT^^' = g^'^Rabcd e\e\ = R^c e\e\ = Ric(e^) ; (18) 

in such a way that the Ricci curvature along any vector of the unitary tangent space coincides with the trace of the 
corresponding 7?. tensor. Incidentally, from the leftmost hand side of eq.([l8|) we see 

14 -\ 

K^c{<,A)^Y^Kf^ ; (19) 

B=0 



that is enough to justify also the last equality in eq. (|1 1]) . Now, given a geodesic passing through q e M, let u e TqM 
its unitary tangent vector. If we put = u, we see that what we called stability tensor for the given geodesic is 
actually 

n\ = , (20) 
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and this completes the proof that TrTi = Ric(u), i.e., = kii{u.). 

As we have just seen, to define the (A/" — 1) sectional curvatures when one direction, u, has been fixed as the tangent 
vector to a given geodesic it suffices to choose the remaining {Af— 1) orthogonal unit vectors {e^}, (A = 1, . . . ,Af— 1), 
transversal to the flow. A preferred set is however obtained if we furthermore introduce the principal sectional 
curvatures for the congruence defined by u, i.e., the local eigenvalues, {A^}, of the stability tensor such that 

n\e/ = X^,e^% (A = 0,...,AA-1) (21) 

with normalized eigenvectors, ||e^|| = 1 ; so, by eq.(^), it turns out 

A„=0 if e„=u. (22) 

and also that the Ricci curvature per degree of freedom kji(u) represents the average of the non trivial principal 
sectio nal curvatures {A^}|^ F rom this last definition, it follows in a perhaps more direct way the well known result. 



e.g.,| Synge and Schild 197^ , according to which the sum of all the sectional curvatures in the {Af — 1) 2-planes 



containing a fixed direction do not depends on the choice of the set of independent normal directions, being simply the 
trace of the tensor 7i^\ In principle, the {Af — 1) non trivial eigenvalues of 7i form a set of indipendent local indicators 
of stability, in what they approximately determine the local behaviour of a deviation vector along the corresponding 



eigendirection; indeed, the analysis carried out for few dimensional systems, [ Cipriani and Di Bari 1997c , has shown 
that the connection between the sectional curvatures and the dynamical behavior do exist, although it emerges 
completely when considered globally, and that estimates based only on naive local analysis can lead only to partial 
answers. It has been claimed that there is a trivial counterexample to this, represented by the fiows on manifolds of 



constant curvature. But it isn't at all a counterexample, as, in that case, Schur's theorem, [ gynge and Schild 1978 1, 
assures that there is no distinction between local and global features of the manifold^. So, although it is generally true 
that the qualitative behaviour of geodesies depends on global rather than local features (and this is the reason of the 
general failures of Toda-like criteria) within the GDA we feel the perception of the intermingled relationships amongst 
them, which allows to gain some global informations from a local analysis of the curvature proprties of the manifold. 

Curvature, Instability and Statistical Behaviour. 

It has repeatedly been claimed that the chaotic and statistical properties of dynamical systems, in the case of mdf , 
e.g. ,| Gurzadyan and Savvidy 198(\ , Kandrup 199C], and for few dimensional systems too, [Szydlowski and Krawiec 1993| 



depend on the scalar curvature of the manifold. Instead, it has been shown with concrete applications, see, for ex- 
ample, [Pettini 1993| or |Cipriani 1993|, that this quantity gives no information at all on the behaviour of trajectories 



for mdf systems, moreover, using the Eisenhart geometrization, it is found to vanish identically, resulting in general a 
very scarcely reliable indicator. The reasons of this failure in both few and many dimensional DS's turn out easily if 
we think over the assumptions implicitly or explicitly made to justify its use. 

The use of scalar curvature has been accounted for neglecting the orientation of the velocity u, so averaging on different 
states of the system, in addition to different orientations of deviation. In [ Cipriani 1993|| the consequences of such an 



approximation are analysed thoroughly; here we show the formal derivation and critically analyse what would be its 
implications. 

If also the tangent vector to the geodesic u is oriented randomly, as ||u|| = 1 too, we can write, in analogy with 
what has been done above for f, except that now we really have Af independent orientations (and then Af — 1 for the 
orthogonal deviation vector): 

K"'> = ^ (23) 

With this strong assumption, we see that the stability tensor loses its dependence on the actual state, forgetting the 
information about the orientation of the geodesic, and the first term in the r.h.s. of eq.(0) becomes 

gbd gac gfj 

" ^^'^'WW^) ^ NiN - 1 ) ^^^^ 



^As we observed, from eq. (g) it follows that the non trivial geodesic deviations should have a component orthogonal to the flow, as a 
parallel perturbation cannot contribute to the possible instability. Stated otherwise, the knowledge a priori of a vanishing eigenvalue of ?i 
suggests to divide the trace by (A/" — 1) instead than M; this is also consistent with the fact that, given a direction, there are only (J^ — 1) 
independent 2-planes containing it. As we were interested in mdf systems, we wouldn't need here to pursue further this point, however 
taking into account such a distinction below, when we average also on the orientations of u, counting correctly the number of (ordered) 
direction pairs as Af{J\f — 1); also because it has to be remarked the great importance held by this factor in the applications of the GDA to 
few dimensional DS's. 

*We will often return in the sequel to some underestimate implications of this theorem, which results of fundamental importance to 
shed light on the properties of the manifold really tied up to the onset of Chaos. 
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where we introduced the scalar curvature of the manifold 

R\ = g'^-Rac = g'^^g^'^Rabcd (25) 
which, using the previous formulae, turns out also to be 

AA-l A/'-lA/'-l A/'-l 

5i.5:Ric(eJ=5:5:'<:-ETr7i'' (26) 
and where we averaged over the right number of directions pahs. Using these approximations, equation (|l^) is replaced 

by 



a z 



ks z^O (27) 

dcf 3? 

where the average scalar curvature has been defined, ks = fcs(q) = //(^J^ — ij' 

We see from the derivation itself that the use of scalar related quantities is unjustified, because the averaging over 
directions of geodesies means to ignore the actual state of the system; i.e., the evolution of the perturbation to a 
given geodesic, along the geodesic itself results to depend on the paths associated to all the geodesies passing on the 
same point q, with arbitrary velocity! This is the same as saying that this claimed average over states turns out to 
be instead an average over 0{J\f^) evolutionary equations amongst which only 0(M) carry meaningful information, 
which turns out obviously suppressed by the 0{M'^ — M) irrelevant ones. 

We will show below that performing that average amounts to assume that the manifold is isotropic; this is definitely 
not true for whatever realistic model of physical DS, whereas it is exactly what happens for abstract geodesic flows 
of mathematical ergodic theory! The only physically meaningful situation in which the scalar curvature carries the 
correct information is that of two dimensional manifolds, i.e., the case in which there is only one sectional curvature, 



equal to the Ricci one, in turn equal to half times the scalar curvature, [Syngc and Schild 1978|. So, any averaging 
process must be carried out carefully, if reliable results are sought. This is true in general, but in particular, for few 
degrees of freedom systems, there are two opposite situations and a list of warnings is mandatory: 



As remarked above, for two dimensional configuration spaces, all the curvatures are directly related each other 
and the results obviously do not depend on the choice made. 

This do not means that for few (but Af > 2) degrees of freedom, the scalar curvature should give reliable 
results. Even more, in these cases, the global constraints in the geometrical features of the manifold causes 
any averaging to destroy the information contained in the sectional curvatures, causing even the Ricci one 
to lead to wrong results, as the averaging can hide all the mechanisms responsible of the onset of Chaos, 



Pi Bari 1996| , [Cipriani and Pi Bari 1997a | 



To highlight the peculiarities of two-dimensional manifolds, we remind that for OS's with two degrees of freedom, 
if interested only in the phenomenological determination of dynamic (in)stability, no matter what geometrization 
is adopted, the results are in agreement. But if the goal is the discovery of the origin of Chaos, it should be taken 



into account that a suitable enlargement of the manifold gives deeper insights, [Cipriani and Di Bari 1997c |. This 



is not necessary if the manifold is already three-dimensional or more, |Di Bari and Cipriani 1997c 



To point out the general (for any J\f) unreliability of scalar related quantities, and even of Ricci curvature in the 



case of small J\f, we recall some results obtained within the Eisenhart geometrization framework, jPettini 1993 1: 
for any conservative DS, 3? = 0, Ric£;(u) = 2 even for a paradigmatic chaotic hamiltonian (the Henon-Heiles 
system). So, eq.(|7|) turns out to be always meaningless, and even eq.(^0|) is completely unreliable for two degrees 
of freedom Hamiltonians. 

Also within the Finsler GDA, it is shown that for a three dimensional manifold (corresponding to a two degrees 



of freedom DS, [Cipriani and Di Bari 1997c[), even the use of Ric(u) leads to incorrect answers. 



The predictions on the N-body system based on the scalar curvature are in complete disagreement with those 
obtained using Ricci curvature, which instead agree with those of standard tools of investigation of chaotic 
properties. 
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The last points relate again to an hasty extension of results borrowed from Ergodic theory of abstract DS's 
and to an inadequate consideration of Schur's theorem. Indeed it has become more and more evident that 
the properties of manifolds associated with realistic physical models differ significantly from those of manifolds 
studied in Ergodic theory and that the instability, when present, has very different sources in the two classes of 
geodesic flows. 

In the latter class, actually, the mechanism of instability resides in the negativity of the sectional curvatures; if 
all the sectional curvatures are always negative, then the flow is exponentially unstable, the system is mixing 



and all the statistical properties are well justified, Katok and Hasselblatt 1995 |. Again, in that case also the 



Ricci and scalar curvatures are always negative, and all approximated/averaged equations predict instability as 
the exact ones do. 

Conversely, the sources of instability for physically meaningful geodesic flows are instead related to the non 
constant and anisotropic features of the manifold. In such a case whatever averaging process risk to destroy 
both non uniform properties, which are moreover each other tightly related via the Schur's theorem. For 
mdf systems simple statistical considerations, actually the use of the Central Limit theorem, allow to recover 
the fluctuations (both in time and directions) and to obtain very accurate answers on the global behaviour, 
[ pipriani 1994 [Casetti and Pettini 1995 , Casetti et al. 1995|. But, when such statistical arguments are not 



justified by low Af values, the use of averaged quantities change completely the answers. 



Incidentally, we remind that the GDA requires some obvious conditions to be applied. We refer to [Cipriani 1993 



Di Bari 1996, Di Bari ct al. 1997, Pi Bari and Cipriani 1997b| for an exhaustive discussion, and recall here the 



apparent singularity of the Jacobi line element, dsj, cq.([^), at the turning points, where U (q) = E. For a 
standard DS, whose kinetic energy is a positive definite quadratic form and then these turning points are only on 
the boundary of the region allowed to motion, it is well known that this do not result in any real singular behavior, 
and, moreover, for mdf systems, the chance that the kinetic energy vanishes becomes smaller and smaller with 
the increase of A/". Nevertheless, the Jacobi geometrodynamics has been applied, pzydlowski and Krawiec 1993 |, 



to a three dimensional General Relativistic DS, for which the kinetic energy is not positive definite. In that case 
the singularity of the metric is inside the region allowed to motion and the Jacobi form of the GDA cannot be 
applied. This is one of the two main points which invalidate the results there obtained, in addition to the use 
of the scalar curvature, whose general lack of reliability (extolled for two dimensional DS's) has already been 
pointed out. It is just in order to overcome this kind of limitations of the Jacobi GDA, that we proposed the 
Finsler extension whose application to Celestial Mechanics and General Relativistic DS's has been performed 



succesfully, [ Di Bari and Cipriani 1997a , Di Bari and Cipriani 1997c [ 



Summarizing, we can say that, while the use of Ricci related quantities is trustworthy for mdf systems, and for 
small Af no averaging procedure is guaranteed to give reliable answers, the scalar curvature has in general, as it is, 
no relationship with the issue of instability, except in the very exceptional situations in which it carries the same 
informations of the Ricci or sectional curvatures. 



Sources of geodesic instability. 

From now on we concentrate on the GDA to Chaos in mdf Hamiltonian systems; therefore we will assume that A/" ^ 1, 

so we can neglect the full set of equations (^ and limit ourselves to the averaged eq. (|l0|) . 

We are left with a formally simple ordinary differential equation for the norm of the perturbation, z{s): 

z" + xis) z = 0, (28) 

where (•)' = d{-)/ds. According to the level of approximation, the quantity x(*) entering in the evolution of z{s) 
would be, in turn: 

Xi-s) = X [q(s),u(s)] = <; (29) 

ksis) = ksHs)] = j^^j^ _ 

From the preceding discussion, it is clear that we will not consider the second case in eq. (p9|) , except when it is 
practically equivalent to the first. Rather, we remind that another set of (A/" — 1) equations which, although being 
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approximated, give more detailed informations on the stability of the flow, could be obtained letting the frequency x{s) 
to assume the value of every non trivial eigenvalue of 7i"c, — (s), with A = (1, . . . , A/"— 1). For systems with few 
degrees of freedom the equations so obtained, yet approximated, have proven to reproduce faithfully the qualitative 
behaviour of solutions of the exact geodesic deviation equations, even when the Ricci curvature fails (nothing to say 



about the scalar curvature), [ Cipriani and Di Bari 1997c |. For high dimensional DS's the problem stands in the huge 



computational work required to diagonalize Ti., rather than in the integration of {J\f — 1) equations instead of one. In 



some cases, however, the stability tensor turns out to be easily diagonalizable, Di Bari and Cipriani 1997c |, and the 
resulting equations help to get some hints on the reliability of the results obtained within the Ricci-averaged approach. 
Equation (^ ) is our starting point for the discussion of the links between geometry and instability. As it stands, it 
points out the importance of the curvature properties of the manifold on the evolution of a generic perturbation to 
a geodesic flow. Later on, we will discuss the steps needed to rephrase the results in terms of the newtonian time t; 
for the time being we limits to general considerations on the stability of the geodesic flow, described as the Jacobi 
s-parameter evolves. 

Most of the results of abstract Ergodic theory refer to geodesic on manifolds of constant negative sectional curvature, 



[ Anosov 1967 , §6.1]. In this case, i.e., when Xa{s) = —a^ = const. , W A = 1, . . . , (A/" — 1), whatever choice is made in 
eq.(|29|), it is apparent that it is also 

X{s) = = const. (30) 

Well then, eq.(|2^) implies that the magnitude of the deviation will increase exponentially fast in s: 

z' 

z{s) — Zo cosh{as) ~\ — 2. gii^}i(Q.5) — > C expajsj . (31) 



In this case it has been rigorously proved, e.g. ,[Katok and Hassclblatt 1995|, that the dynamics possess the strongest 



statistical properties, as positive Kolmogorov entropy, mixing, exponential decay of correlations, and so on; and it is 
in fact the most chaotic one. Therefore, if a DS were characterized by constant negative sectional curvatures, then it 
would be strongly chaotic and approaching the equilibrium state within a relaxation time of order Sr ~ Q;^^, or, in 
terms of the newtonian time^, 

Tr « i{W)t ■ (32) 

As far as we know, in the case of non constant but everywhere negative sectional curvatures, 

Aa(s) < -/3' < 0, VA=1,...,(AA-1); Aq = , (33) 



while the flow is deflnitely unstable, |]Katok and Hassclblatt 1995 1, for what concerns the statistical properties, it is 



plausible, and there are convincing argumentations, according to which also the approach to equilibrium take place 
on a s-time scale of order /3~^, although we are not aware of any rigourous and explicit theorem. 
Bearing in mind that the exact equation for the evolution of geodesic deviation is eq.(|^), of which eq.(p8|) is an 



approximation, we should discuss at this point the relevance of Schur's theorem (see, e.g., |Synge and Schild 1978 
§4.1]), which refers to isotropic Riemannian manifoldsp^t whose deflnition is the following: 

"An A/'-dimensional {Af > 2) riemannian manifold, M^, is said isotropic at a point q e if the 
sectional curvatures if'^^^(x, y) evaluated at q, do not depend on the pair (x, y), i.e., on the 2-plane." 

In this case the Riemann tensor has the form 

Rabcd = IC {gacgtd - gadgtc) (34) 

where /C is a scalar quantity, a priori depending on the position, /C — /C(q). Well, then the Schur's theorem assures 
that: 



more careful evaluation of the physical relaxation time would give the result 



<a-Tr< {W)y 



Although the two limits do not differ significantly for standard mdf systems, as the FPU chain, in the case of very inhomogeneus gravitational 
N-body system out of global virial equilibrium, the difference can be quite remarkable. This has no consequences in this context because 
both them have neither constant or negative sectional curvatures. The implications of these mo re accurate limits will be discussed elsewhere. 
'^''The extension to more general differential manifolds, as the Finsler ones, is discussed in |]Rund 1959 . 
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"If an A/'-dimensional (J\f > 2) riemannian manifold, Af/v', is isotropic in a region, then the riemannian 
curvature is constant throughout that region" 

That is, the scalar /C is a constant quantity, do not depends on q. This incidentally implies also that the covariant 
derivative of Riemann tensor vanishes. 

What this theorem entails on the issues we are facing on can be deduced looking rather at its reverse implications. 
Indeed, according to the theorem, if a manifold is isotropic, then all the sectional curvatures are also constant; 
nevertheless, if the manifold is anisotropic, nothing can be said about the constancy or variability of the sectional 
curvatures, as it could be eitheip] 

ii:*^^)(xi,yi) = ai = const. 7^ if^^^ (xa, 72) = 012 = const. , 

or ^^(^^(x, y) = A{s), a variable quantity; and in this case the theorem do not give more informations. Then, this 
occurrence can be interpreted saying either that Xa 7^ in consistent with both variable or constant {A^i}, or, 
conversely, that Xa = const., V A is consistent with both isotropy or anisotropy. 

Other interesting back implications of the theorem explain what can be learned looking at the behaviour of Ricci 
curvature. 

• If at least one of sectional curvatures varies, Xa — Xa{s), then the manifold cannot be isotropic. 

• For an isotropic manifold, i.e. for Xa = a, VA — 1, . . . , {Af — 1), then also the Ricci and scalar curvatures are 
also constant, and — kg — a. This is the sole situation in which the averaged quantities equal exactly the 
the sectional curvatures. 

• If the manifold is anisotropic, then the Ricci curvature along different geodesies assumes distinct values, whereas 
the scalar curvature do not distinguish between various initial conditions. 

• Even in an anisotropic manifold, if the sectionals are constant, then both kji and ks do not vary. In this case 
they represents their averaged values. 

• When at least one of the sectionals varies, the manifold must be anisotropic, and in this case the averaging 
procedure can hide the fluctuations: with fluctuating {Aa(s)}, it is possible to have either kii{s) also vary- 
ing with s, or even constant. These remarks apply with even more relevance for scalar related quantities, 
[|Cerruti-Sola and Pettini 1995 . 



• In particular, if we find that fc/j = const., we cannot say nothing about the behaviour of sectional curvatures, 
because this is consistent with all the possible features, namely, we can have either an isotropic or an anisotropic 
manifold with constant sectionals, or even an anisotropic one with fluctuating Ay^'s! 

• Vice versa, if the outcome is that kii{s) is a fluctuating quantity, we can definitely assert that also the sectional 
curvatures vary along the geodesic, and, by the Schur's theorem, that the manifold is actually anisotropic. 

Luckily enough, almost all physically meaningful geodesic flows fall into this last category, and either farther from 
complete integrability we are or higher is A/", more true it is. According to this picture, the manifold is anisotropic, 
with sectional curvatures rapidly fluctuating, almost independently from each other, in such a way that, when M 3> 
1, we can assume (and even verify!) that Ric(u) is the sum of A/" — 1 uncorrelated quantities, of which kfi{u) 
represents the average. For mdf systems, and when the motion is definitely far from quasi-periodicity, it can be 
shown, [Cipriani 1993, Casetti and Pettini 1995, [Casetti et al. 1995 |, that it is consistent also to assume the random 



character of the quantities entering Ric(u), and using the central limit theorem, to relate the fluctuations of kf;{u) 
to those of the Ay^'s, obtaining an analytical estimate of the asymptotic rate of growth of the norm of the deviation 
vector, i.e., of the LCN. 

Once realized that the case in which all the sectional curvatures are constant is an exceptional one, and consequently 
that the rigorous results on abstract geodesic flows carry very little physical interest, we turn to the study of the 
evolution of a perturbation to a realistic geodesic flow. 

^^In the sequel we will rephrase these argumentations in terms of the principal sectional curvatures along the flow, (the Ayi's), though 
these depend on a parallely propagated vecto r. This is legitimate as long as H satisfies very generic hypotheses. A very interesting account 



on these general arguments can be found in [ Berndt and Vanhcckc 1992 1 
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Mechanisms responsible of the onset of Chaos. 



Bearing in mind the warnings made above, we now look at eqs.(|2^) and (|2S|), in order to see under what conditions 
instability can arise. It is obvious that when — —o? — const., these equations predict an exponential growth of 
the magnitude of the geodesic deviation. As stressed before, a constant value of Ric(u) can occur for very different 
behaviour of sectional curvatures, but a negative constant value implies that the average value of them is permanently 
negative, and this do not happen for any geodesic flow of physical significance, neither it found that fc^ is mostly 
negative. To see this, let us write the expressions for Ricci and scalar curvature in the Jacobi metric for a conservative 
DS with M degrees of freedom, [ pipriani 1993 : 



and 



A:5(u) 



NU 



(gradW)^ 



ds 



d%l_ 

ds^ 



(35) 



(36) 



where we indicated the usual Af-eucUdean Laplacian and gradient operators with A and grad, respectively. A first 
look at the expressions shows that, on passing from the average over directions of the deviation vector to the loose 
approximation given by the scalar curvature, some important information get lost. Indeed in both kp. and ks appear 
the A/'-dimensional laplacian of the interaction potential, with a weight reduced by a factor 2 in the former, and the 
A/'-dimensional squared gradient, i.e., the sum of the squared forces acting on each particle of the system. One main 
difference is on the sign and on the weight of this term, which is amplified by a factor 0{N) in the scalar related 
frequency, ks to which it contributes negatively, thus increasing the chance of trivial instability of the solution. At 
least equally relevant to the evolution of the instability is the lacking in ks of the last two terms in fc^j, as they 
describe the explicit s-time dependence of the curvature of the manifold, in particular for systems out of global virial 
equilibrium, i.e., subject to collective oscillations which involve macroscopic fluctuations of total potential and kinetic 
energies, as it occur for a collisionless N-body system during its very early stages. The consequences of the differences 
just listed can be grasped intuitively if rephrased in terms of the newtonian t-time. As a matter of fact, we see that 



the last two terms in the r.h.s. of eq.(35) are nothing else than the logarithmic time derivatives of the kinetic (or 
potential) energy: 



1 



dU _ 
'dl ~ y/2E-U 



1 W 



ds^ 



1 d W 
'2Wdi\W 
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2VV 




(37) 



and, for a mdf system, these relative fluctuations arc damped by a statistical factor that obviously increases with Af. 
Moreover, they exploit an explicit dependence of fc/j(s) on rapidity with which the global virial equilibrium is attained. 
The magnitude of those relative fluctuations indeed, though generally smaller and smaller as the number of degrees of 
freedom grows, depends in general on the virial ratio, and could be not negligible, in particular for a self-gravitating 
system, where collective effects are likely to occur, due to the long range nature of the interaction. Thus, whereas 
the Ricci curvature keeps memory of these evolutionary processes, the scalar one, just because forget ab initio the 
real dynamics of the flow, erasing the u dependence, cannot take care of the consequences of the approach to this 
dynamical global equilibrium. If we add to these remarks the presence of (Af — 1) Ricci curvatures along directions 
extraneous to the dynamics, the wrong sign and weight it attributes to the square average forcepl, with respect to 



Ricci curvature in both Jacobi and Finsler cases, Cipriani 1995 , Pi Bari 1996 , Pi Bari et al. 1997 1, together with the 
evidence of the effectiveness of the use of the latter in gaining deep insights on the sources of instability and in the 



analytical computation of instability exponents, [ Casetti et al. 1995 , Cipriani 1993 , Cipriani and Di Bari 1997c |, it 
turns out once more why the use of scalar related quantities is ever more unreliable in the case of mdf systems. 
A first look to the expression of the average curvature in the {Af — 1) 2-planes containing u leads to the following 
considerations: 

• The A/'-dimensional laplacian is almost everywhere positive for every confining interaction. It is always positive 
near the minimum of any potential and can be locally negative only in limited regions. 

• The squared A/"- gradient contributes with a positive sign to k^. 



^■^Which has been responsible of a largely unjustified debate in the past, |Gurzadyan and Savvidy 1986, Kandrup 1990[ . 
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• For J\f > 2 also the square of the of the s-time derivative of the potential energy contributes with a plus sign. 



Then, the term with the second derivative of U alone can give a negative contribution to the Ricci curvature. 



What it has been found for the FPU chain, either numerically integrating the trajectories, [ Pettini 1993| and jCipriani 1993| 
or canonically averaging over phase space, [ Casctti and Pettini 1995| is that the Ricci curvature (and then kji) is 
mostly positive, that its average value is always positive, that there is no correlation at all with the degree of 
stochasticity in the dynamics (as measured by, e.g., the maximal LCN), |Cipriani et al. 1996 |, that the chances to 
find a negative value for it decrease quickly increasing Af, vanishing in the thermodynamic limit. Moreover, in 



[Cipriani 1992], it has been found the numerical evidence of the actual occurrence of the foreseen virial transition 
process, [ Cipriani and Pucacco 1994a , Cipriani and Pucacco 1994b[ . In the figures |l| and ^, this transition is neatly 
detected, as well as its dependence on the number of degrees of freedom and on the regime of chaoticity. It is found 
that, while during the first phase of virialization, the frequency of events such that < is not negligible, after then, 
and already for relatively small values of A/" 10^), the probability p of a negative value for ku, is actually vanishing. 







p{kR < 0) - C^fcxp[-{t/T,f]-^^^^ 

as highlighted, from the pure behaviour of the measured cumulative frequency, _F_ {t)f^. 
dcf N^it) _ iV(/cfl <0) 



(38) 



(39) 



N{t) N,teps 

shown in figure |l|. We observe, moreover the strong decrease of this frequency with Af, and the very weak correlation 
with the degree of stochasticity, as shown by figure |^. The following discussion explain the qualitative behaviours of 
the constant CV and the virialization time Ty = Ty{Af, f3e). When it was firstly detected, it was surprising enough to 
see the lack of correlation between the degree of stochasticity, as measured by the maximal LCN, or the parameter (3e, 
and the frequency of occurrence of negative values. Even the Af dependence was contrary to what expected on the 
basis of the belief that the statistical description is more reliable as the number of particles increases. On the light of 
the analytical estimates, since then these results turn out however completely understandable. A thorough analysis of 



these results is presented in [Cipriani 1993 1 and completed in [Cipriani and Di Bari 1997b , but we can say here that 
the efficiency of the virialization process is clearly decreased by the occurrence of a quasi periodic behaviour, i.e., by a 
nearly integrable dynamics. Moreover, the phase mixing, which is in turn responsible of the trivial loss of correlations 
between particles, becomes faster along with the increase of the largest normal mode frequency excited, i.e., with Af 
too. 

As the unpredictability is believed to growth with the numbers of interacting particles, the inverse correlations between 
Af and p(fci? < 0), once more force us to discard the hypotesis of an instability driven by negative values of Ricci 
curvature. Besides, the results obtained within the Jacobi picture are in complete agreement with those coming out 



from the Eiscnhart geometrization, [ Casetti and Pettini 1995 , Cipriani 1993 |; and in this last framework, the Ricci 
curvature for the FPU chain is always positive (more, it is always kj^-^ > 2), nevertheless, is there possible to recover 
in a very elegant and effective way all the (in)stability properties of the dynamics, obtaining an analytical algorithm 
for computation of the largest LCN. 

If the sign of k^ is mostly positive, and seeing that numerical integrations of eq. (p8|) allow to recover all the qualitative 
and quantitative stability properties of mdf systems, we are led to ask where stems from the exponential growth of 



the geodesic deviation. The answer is actually written in any textbook, e.g., [Landau and Lifsits 198C[, and resides 
in the mechanism of the swing: the instability is driven by parametric resonance induced by the fluctuations of the 
positive value of the frequency x(s). We refer to [Cipriani 1994 [Pettini 199^ [Casetti and Pettini 1995| , [Di Bar i 1996 
Di Bari et al. 1997 [ for details, and here simply recall that in the non-autonomous equation 



■ + n^{t)x = 



(40) 



with 



n'{t) = nl[i + f{t)] 



instability (i.e., exponential growth) of the solution can arise, fixed the average frequency Oq, if suitable conditions 
are fulfilled by the modulation factor f{t). In the case of eq. (p8[) , the modulating factor is not an explicit function of 
s-time, rather depends on s through the actual state [q(s), u(s)], of the system. So, in order to foresee the behaviour 

Defined as the ratio between the total number of occurrence of negative values of Ricci curvature and the total number of steps 
performed. That is, A^stcps = t/At, At being the integration time step. 
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of the deviation vector it is important to know not only the average value Xo '= (x('S))s, but also the amplitude, 
'^x ~ ((x(*) ~ XoY)s and the variability time-scale Ts of its fluctuations. 

Once again, we see why the cancellation of the fluctuating terms consequent to the averaging over u, is at the origin 
of the unreliable results obtained with scalar curvature. 

In order to compare the outcomes obtained within the GDA with those resulting from standard dynamical system 
approach, we need to rewrite the equations in terms of the newtonian i-time, exploiting the non afflne relation given 
by eq.(Q]). If we do this, equation (BS) becomes: 



(41) 



dcf 



where kn '== 2yV^ kn is a sort of rescaled Ricci curvature. We are left wiht an ordinary differential equation for 
a damped harmonic oscillator, whose peculiarities lie both in the fast variability of the frequency kji{t), and in the 
undefined sign of the damping term, which actually has zero average and can act alternatively also as an anti-damping. 
The presence of this term is nevertheless crucial in order to understand the sources of instability in the pliysical time 
gauge. Indeed, an equation completely equivalent to ( p| ) is obtained with the change of variable 



Y{t) ""^^ W-i/2(i) ^(i) ^ 



(42) 



where it is understood that all the quantities depends on t through the one-to-one correspondence s = s{t), fixed by 
ds = ^/2yV dt, yV > 0. For mdf systems indeed the probability that W = is negligible, and, for a confined DS, the 
virial theorem guarantees also that W{t) is a quasi periodic function, bounded from abovepi; in such a way that the 
qualitative long term time behaviours of Y{t) and z{t) are the same. This well known substitution allows to get rid of 
the damping term, leading to a standard Hill's equation: 



f(t) + Q{t) Y{t) = , 
where the frequency Q{t) follows from the rescaled average curvature kji{t): 



Q{t) 



Q{q(t)} knit) 




1 W 

2 W" 



(43) 



(44) 



As heralded above, the terms appearing in the expression of Q{t) with respect to knit) are very important for the 
understanding of the behaviour of the perturbations in the physical time. To see why, let us write down the explicit 
f-expressions: 



knit) = ^ 



(gradW)^ 



+ (AA - 2) 




1 W 
2W 



(45) 




(46) 



We see that the terms arising from the s derivatives of the potential energy, and here rewritten in terms of i derivatives, 
appear in Q{t) with a weight that, in the large N limit, is actually reduced by a factor JV^^. The relevance of this 
outcome becomes self-evident as we consider the the sources of instability in eq. (^ , or (|4l|). Indeed, as long as 
we consider customary mdf systems, as emphasized above, it is never found that Q(t) < 0, and the frequency with 
which knit) < 0, is always very small, becoming absolutely negligible with increasing Af, in which case asimptotically 
vanishes, as the virial equilibrium is attained; see figures |l| and ^ and also |Cipriani 1993, §5.2.1]. In addition, this 
frequency do not show any correlation with the degree of stochasticity, as it is evident from figure ^, obtained from 
[Cipriani 1992], see also [Cipriani et al. 1996]. 

The occurrence of Chaos, i.e., the exponential growth of all the solutions (except, possibly, for a set of initial conditions 
of zero measure) of eq.(^3|) is consistent with a positive sign-definite Q{t) if the mechanism responsible of the onset of 
instability is the parametric resonance mentioned above. It is crucial for this phenomenon to occur that the amplitude 

'^^This is not strictly true asimptotically for the mathematical gravitational N-body system, but it is true for almost all initial conditions 
at any finite time. Moreover, and this is an important point, W is always bounded from above, when the physical interaction is addressed. 
We stress that the strategy we introduce below to cope with the mathematical singularities of the gravitational potential do not change 
its very nature of non-compact DS. 
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of fluctuations fulfil some conditions, whose details depend either on their quasi periodicity or on their stochastic 
behaviour. This explain the relevance of the reduced weights in front to the fluctuating terms in Q{t) with respect to 
kii{t) in order to recover the known results in terms of the newtonian time. 

The argument just discussed enlighten also why a mdf system should go across a transition in its qualitative behaviour 
along the approach to the global virial equilibrium, during which the amplitude of the fluctuations in its macroscopic 
parameters, as the kinetic and potential energies is damped just in consequence of the trivial phase mixing previously 
recalled]^. When the virialization phase is accomplished, the instability of the dynamics is governed mostly by the 
fluctuations in the dominant term, which, accordingly to a kind of principle of equivalence of any geometrization (in 
any case confirmed by analytical calculations and numerical computations) is the A/'-dimensional euclidean laplacian, 
see table 0. 

Being the translation of the dynamics in (Jacobi) geometrical language almost straightforward for any additive 
interaction potential, U, we could achieve order of magnitude estimates of all these quantities in the general case; 
nevertheless, to avoid a (still) heavier discussion, we will consider separately the FPU unidimensional chain and 
the gravitational N-body problerrf^ that can be considered most representative among the short and long range, 
respectively, interaction potentials. 



Standard potentials: the FPU chain. 

The first mdf system investigated by computer simulations has been the celebrated unidimensional chain of N weakly 
anharmonic oscillators. According to the power-law of the anharmonic terms, cubic or quartic, these models were 



indicated as FPU-a and FPU-/5 models, respectively, [ Fermi ct al. 1965 1. We will focus our attention here on the second 
kind, essentially because it has been thoroughly investigated, both from dynamical as well Statistical Mechanical points 
of view, see, e.g., [Pettini and Cerruti-Sola 1991, Benettin et al. 1984| and references therein. 
In this case d = 1, so Af ^ N and the Hamiltonian is 



N 



^(x,p)-E 



1 2^1^ 



(47) 



We recall that, in the quasi- harmonic limit, namely for /3e <C 1, where e E/N is the specific energy {E = _ff(x, p)), 
the introduction of the normal modes allows to exploit the properties of this quasi-integrable system. Such a good set 



of coordinates, (X,P), is defined as, [ Toda et al. 1992 1: 

Pk = Xk 



Xi = \ — 



1/2 JV-l 

E 

fc=i 



Xk sin 



(48) 



and is well suitable for a perturbative treatment. The Hamiltonian, eq.(47), in the new coordinates reads 



N 



k=l 



N 

E 

{jl J2 J3} = 1 



C{k,juj2,j3)XkX,,X,,X^ 



33 



(49) 



with ujk = 2sin(-^) and the C(fc, ji, j2j js) are coefficients depending on the choice of the boundary conditions (fixed 
or periodic). Using these coordinates the harmonic limit, /3e — *■ 0, apparently turns out to be integrable (actually N 
independent harmonic oscillators!), and the presence of the anharmonic term is responsible for the coupling of the 
normal modes, and the loss of integrability. 

However, since we are interested in the behaviour when the system is very far from integrability, the use of the normal 
modes coordinates, although enlightening for the analytical estimates at low e, is of no help in the strong stochasticity 

By this expression we mean the purely kinematic effect of particle-particle correlation loss, entailed by the different orbital times, 
not to be confused with the mixing property of Ergodic theory. This explain also because very near to integrability, being the motion 
quasi periodic, such correlations persist in time, decreasing even the speed of this phase mixing, which is in turn responsible of the 
convergence of the curvatures to their asymptotic values. We already discussed the increasing effectiveness as the thermodynamic limit 
is approached. These remarks on the relevance of Af are intended to warn about the risks hidden in extrapolating naively the outcomes 
of numerical simulations performed with too small Af; in particular for non extensive interaction potentials, like the gravity, where, as we 
will show, the A/'-dependence is often hardly predictable. In these situations, the numerical evaluation of average values and dispersion of 
geometrodynamical quantities should be performed with long enough runs involving a suitably extended J\f range, to single out possible 
scaling behaviour. 

^^We will restrict to the physical case of three spatial dimensions, although many theoretical hints on the issue of relaxation in N-body 
systems have been achieved studying two or (mos tly^ one-dimensional systems, which display the very nice property to allow for computer 
simulations virtually free from numerical errors, [Fsuchiya et al. 1996 1. 
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regime, and it is superfluous from a numerical viewpoint, as it would imply a reduced efficiency of the algorithm. 
To estimate the order of magnitude of various terms entering the Ricci curvature, we start inquiring on the dynamical 
time scale of the system. This is relevant to our analysis as it enters in the s derivatives of the potential or kinetic 
energies. As most of interaction potentials, the FPUt3 model possess in the genoma its own lenght and time scales. 
These generally depend on some global parameters describing the macroscopic state of the system (as temperature and 
density). This DS is however peculiar in this respect being isochronous at low energy, as any harmonic system; i.e., 
when Pe <^ 1, the global dynamical time-scale is of order unity. An important and well known property of this model, is 
the existence of a strong stochasticity threshold (SST) as the specific energy e increases above a critical value Cc, being 
(3 held fixed, distinguishing among two different regimes of Chaos, and also involving a transition on the features of the 



relaxation processes driving the system towards the equilibrium state, [ Pettini and Cerruti-Sola 1991 , Pettini 1993 |. 
This SST is easily located by the intersection of two different scaling law of the maximal LCN, 71, with e: the first 
results, reported in the references above, seemed to support the claim that 



71(e) cx 



for 



71(e) cxe^/^, for 



(50) 



where ec ~ A series of numerical computations with a large enough N and carried out for very long integration 
intervals, |Cipriani 199S |, extended the previously investigated range of specific energies, and while confirming the 
power law at low energies, find out that above the SST, the scaling is 



71(e) cx e^^^, for e >, 



(51) 



These outcomes are in complete agreement with the elegant SM canonical calculations, performed in the thermody- 



namic limit, in |Casetti and Pettini 1995, Casetti et al. 1995 



Within this context the GDA revealed is effectiveness in pointing out this threshold in a straightforward manner: figure 
^ shows the behaviour or (fcj^(e))s for different values of anharmonicity. The transition energy is neatly detected and 
it should be remarked the definitely quick convergence of the Ricci curvature to the average values there reported, for 
all values of f3e, as it is shown in |]Cipriani 1993 , §5.3]. 

Moreover, within the GDA, it has become feasible to use a semi-analytic computation of the largest LCN, performed 



in [ Cipriani 1993 | using the computed time average, dispersion and correlation time of the Ricci curvature, and in 
[ Casetti and Pettini 1995 , using the same quantities, evaluated instead from a SM point of view, averaging over phase 
space. 

The corrected scaling law behaviour found above the SST, though perhaps not crucial for what concerns the exis- 



tence of two different regimes of Chaos, is relevant to our present discussion in that a simple calculation, [Cipriani 1993|, 
shows that in the strongly anharmonic regime, the dynamical time scale, to, of the FPU-/3 model scales exactly as 
e""'^/^, in such a way that the quantity Li '= 71 • tu is constant when strong Chaos has developed. This point is of 
outmost importance for the understanding of the nature of instability present in the gravitational N-body system. 



In [ Cipriani 1993 [, where it has been explored the dynamical behaviour of the FPU-/3 system also for a range of values of 
the anarmonicity parameter, /3, it has been shown that the relevant quantity that determines the level of stochasticity 
of this DS is, as expected, the product /3e. Therefore, in the sequel, we will assume that the anharmonicity parameter 
is fixed (namely, /3 = 0.1) and let e to vary. Moreover, all the quoted numerical results refer to simulations performed 
adopting periodic boundary conditions, i.e., xq = x^, x^f+i — xi. 

Once established the scaling behaviour of the dynamical time with energy, we turn to the order of magnitude 
estimate of the quantities entering the Ricci curvature for the quartic FPU chain. 
For example, the explicit computation of the A/'-laplacian is straightforward and yields^ 

2 + Tf - ^»)' = 2 + ^mxf)u , (52) 



AU 



which can be also written as (for the details, see [ Cipriani 1993 , Cipriani and Di Bari 1997b| ) 

AH r / / 

1 + 3 Vl + /3e-l 



to 



whereas for what concerns the dynamical time scale it is easy to find that 

ni/2 

X(c) 



Vl + 4/3e- 1 



2/3e 



(53) 



(54) 



Incidentally, this qu antity also coincides with the Ricci curyature per degree of freedom in the Eisenhart metric, cfr. 



Casetti and Pettini 1995 
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where 1(c) is almost constant dimensionless integral, and c is a parameter measuring the departure from harmonic 
behaviour, and so is directly related to fie: c = 1 in the harmonic case, c = in the purely quartic domain. It has to 
be remarked that 1(c) is a monotonically increasing function of its argument, whose extremal values are 

1(0)^— ^ =1.311... ; X(l) = - 1.5708... 

4[r(|)]^ 2 

Similar expression can be found about the (A/", e) dependence of all other quantities, as resumed in table |^. In 



[Cipriani and Di Bari 1997b|, where a comparative analysis between FPU-/3 model and self gravitating N-body system 
is also presented, we propose an explanation of the scaling law behaviours of Lyapunov exponents on the basis of these 
estimates, supplemented with analogous evaluations on the amplitude of fluctuations, during and after the virialization 
process. 

A preliminary remark is devoted to the relation of quantities entering the frequency Q{t) (or kr{t)) to the dynamical 
time to- While at intermediate e values the relationship can't be singled out by eyes, in both the asymptotic regimes 
it is easy to see, that all the quantities behave as tj^^, the different weight they have into Q{t) being solely due to 
their Af dependences. We see from table ^that among the static terms^, the one related to the average squared force 
is always depressed by a factor A/""" with respect to the laplacian, with at least a > 1. The use of the central limit 
theorem allows then to state that, for the FPU-/? chain, the fluctuating terms (see note ^8|), although comparable 
with the gradient term when the system is out of virial equilibrium, are then always negligible, in the large N, with 
respect to MA/M, and are even damped by a further factor of at least 0{Af~^^'^) or 0{Af~^) when that equilibrium is 
attained. As the term for which the damping is less effective, i.e., that with the second order derivative with respect 
to time, has also a zero time-average^ we understand why this it is important to fully appreciate the relevance of this 
transition in mdf systems. 

On the grounds of the above discussion, we understand why the results obtained within the frameworks of the 
Jacobi and Eisenhart metrics agree for the FPU-/3 model: the geodesic instability is driven by the fluctuations of the 
GDI's rather than by their negative values. For the FPU chain, all the terms which appear in Q{t) in addition to 
AU/Af, which is nothing but the Ricci curvature per degree of freedom in the Eisenhart geometrodynamics, after the 
virialization phase has been completed can be neglected not only with respect to the average value of the latter, but 
even with respect to its own fluctuations, as the geodesies explore the manifold. This is generally true at all regimes. 



[Cipriani and Di Bari 1997b |, and is always true above the SST, where the analytic computations of the maximal 



LCN obtained using some elementary results of the theory of Stochastic differential equations, [Van Kampen 1976|, 
coincide, irrespective of the metric used. For small values of either f3e or A/", the quasi-periodic nature of the motions 
as well as the reduced effectiveness of the damping cause a much slower convergence to asymptotic values of the 
fluctuations, so reducing the reliability of flnite time estimates. 



Long range, unscreened interactions: the gravitational N-body system. 

When we focus our attention towards the gravitational N-body system, we are at once faced with an ambiguity 
which cannot be solved if considered only from a mathematical point of view. Along with physically relevant reai 
peculiarity indeed, the Newtonian Gravity possess some mathematical singularities which do not actually involve any 
physical singular behaviour. So, in translating the dynamics into a geometric picture, we will try to cope with these 
mathematical divergencies without modifying the physically peculiar features of this interaction. 
To describe the dynamics of N bodies of mass {mi} interacting via the gravitational potential 

^(^ii) = -Gmimj/rij , 

we introduce the coordinates {q"} in the 3N-dimensional manifold AI^ 

M : {g"} t.c. U{q) < E{< 0) ; 

where: 

= ^/^rf ; a = 3{i ~ 1) + a and i = I, . . . , N; a = 1, 2, 3 =^ a = 1, . . . , 3iV = TV. (55) 

^®This is only a rather generic terminology, to let it be understood that these terms do not involve explicit time derivatives. Of course, 
they also slowly evolve along with the process of approaching the global dynamical equilibrium, though in a very different fashion in the 
two DS's we are considering. Nevertheless, the virial transition is much less evident for them than for the other two terms, in such a way 
that, in the FPU-/3 chain as well as in the gravitational N-body system too (although with a less clear hierarchy), after the first transients, 
the fluctuations responsible for the onset of Chaos come essentially from them alone. 
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endowed with the Jacobi metric, eq.(|^), with the total potential energy given by 

N-l N N 

^(q) = E E ^(^^'») ^ 9 E E ^(^'^) (56) 



2 



In order to get rid of unnecessary complications, we assume that all the particles have unit mass^^ in such a way 
that the total mass of the system is N, its average density p = N/V, where V is the volume occupied by the points. 
Amongst different possibilities, we selected this choice of units because it seems to be one of the most reasonable, and 



moreover allows for a direct interpretation of the results recently obtained, | Cerruti-Sola and Pettini 1995 1. 
In analogy with what has been done for the FPUtS model, we start to inquire about the possible existence of any 
scaling law of the dynamical time, tjj, for the N-body problem. This question is a very crucial one, because of the well 
known scale-free nature of Newtonian Gravity. Nevertheless, an N-body system do possess indeed a natural time-scale, 
which is known in stellar dynamics as the orbital time, because it represents a suitable average of the typical period of 
a body moving in the system. Neglecting some numerical factors of order unity, completely irrelevant for the present 
discussion, this dynamical time-scale is usually quoted in terms of the average mass density: 

tn oc i?3/2/^i/2 (57) 

where D ~ V^l"^ represent the typical spatial extension of the system, and we have exploited our choice on the mass 
spectrum, neglecting again constant factors. It is easy to show that tn represents the right order of magnitude for 
most orbital periods even in a moderately inhomogeneus system, except for the small fraction of tightly bound core 
objects. 



In a recent paper, [ Cerruti-Sola and Pettini 1995 1, where for the first time were performed numerical simulations 
specifically devoted to support (or even correct) the GDA to Chaos for self-gravitating N-body systems, it has been 
reported, along with other interesting analyses, the absence of any threshold in the behaviour of the maximal Lyapunov 
exponent along with the increase of the energy density e '= E/N. Moreover, it has been there discussed the seemingly 
non physical behaviour of this system when the number of degrees of freedom is varied. It has been found, indeed, that 
increasing the number of particles, the degree of stochasticity, as measured by the largest LCN^, decreases, contrary 
to physical intuition. Statistical Mechanics expectations, and to what observed in any other mdf system considered 
before. The authors refused to believe in this apparent discrepancy, and attributed the cause to the too strong 
approximations which, starting from exact JLC equations of geodesic spread, led to an equation like our eq. (^3|) . We 
now show that their numerical results are indeed correct, and that no failure can be attributed to the averaging process 
that lead firstly to equation and down there to eqs.(p8|) and (|4^). The seemingly unphysical result, turns out to be 
instead completely reasonable, if we consider how the numerical simulations were performed]^ As remarked above, 
all the bodies have been assigned unit mass, and initial conditions are generated from a spatially uniform distribution 
whose lenght scale D is chosen to give the desired value of the total energy E, being the velocity (components) of the 
masses populated according to a gaussian distribution (in the CM coordinates system) whose variance, tr^, is adjusted 
to fulfil the virial theorem, i.e., 2 (T) — so that D and N completely determine the value of total energy, 

E = {U)/2 — >U/2. 

With this choice of initial configurations, in the case N 1, the total energy {G = = 1), 
7V ^ N-i N 

i=l k=i+l 

scales obviously as 

where Fp is a factor of order unity depending on the actual density distribution, and it is so allowed to evolve slowly 
with time, along with the development of a typical core-halo structure, although this happens on time scales longer 

^^The problems we are now discussing do not depend on this point. It is obvious that this unrealistic assumption should be removed 
when issues as the problem of equipartition of energy, and the consequent mass segregation is addressed. But the conceptual question we 
are focusing on now would only be made heavier if a non uniform spectrum of masses is allowed. 

^"^For a non compact phase space, as that of the n pwtonia.n N point masses, it is not rigorous to speak about true LCN's, but we can 



equally adopt this terminology, referring to the paper, | Ccrruti-Sola and Pettini 199f| for a critical discussion of this item. Besides, we defer 
to the following discussion and to a forthcoming paperT ^ipriani and Di Bari 1997h[ |, where some of the open questions will be addressed 
in a hopefully self consistent way. 

^^Actually these scaling laws, reported in table ^, do not depend qualitatively on this particular choice. 
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than those lasted usuaUy by numerical simulations. So the total energy scales essentially as \E\ w N'^/D. What this 
implies on the relationship between specific energy e and the dynamical time scale t^i? The answer is straightforward: 
we have |e| w N/D; on the other end we have also 

_ py^ {Ny\E\f/^ _ N 

- ^1/2 ~ - |g|3/2 ' ^""^ 

which is a result valid at any energy. Incidentally, we note that actually Id ^ D/ay. This simple calculation explain 
all the results found in | Cerruti-Sola and Pettini 1995|; indeed the scaling behaviour of the largest LCN as |ep^^, with 



held fixed, is consistent again, as for the FPU chain above the SST, with a constant value of the quantity Li = to- 
Even the (un)believed unphysical result of a decreased stochasticity when N increases, turns out to be consequent 
to the neglecting of the N dependence of tu, and its interpretation becomes instead very clear within this picture: 
again, the dynamical time is proportional to N, if e is held constant, in such a way that the instability exponent 
decreases in absolute value, but remains constant if measured in units of the dynamical frequency. In a sense, the 
lack of stability property, [ Ruelle 1988| , by the Newtonian potential, i.e., the non-extensivity of the potential energy, 



reflects itself on a dependence on N of the dynamical time scale which do not occur in any customary interaction. 
The preceding discussion so contributes to give to the results presented in Cerruti-Sola and Pettini 1995|| a physical 



interpretation which largely confirms what argumented there and moreover support in a stronger way the use of the 
GDA, which turns out to be very fruitful even in the application to this peculiar and singular mdf system. It has been 
removed in fact the doubt about the reliability of the averaged equation (^, which really carries out almost all the 
informations about the qualitative behaviour of the dynamics of mdf Hamiltonian systems. The scaling law behaviours 
found numerically agree with a surprising precision with the analysis here proposed, for the LCN computation, and 
also for the evaluation of some quantities, named Geometric Chaoticity Indicators {GCI) by the authors: we suggest 
to look at figures 4 and 6 of the cited paper, to see how both the e and N dependences are perfectly fitted by the 
analytical estimates here presented (we remark that the scales used in the figures there are in natural logarithms). 

Nevertheless, when we try to apply the GDA to gravitational N-body system, another issue should be necessarily 
addressed: whatever geometrization procedure is adopted, the term in the Ricci curvature (or even in the scalar 
one) which is usually dominating over the others has in this case a very singular behaviour. The laplacian of the 
Newtonian potential indeed vanishes everywhere except on the positions occupied by the sources of the field, where 
it diverges. As one of the most successful applications of the GDA for high dimensional DS's resides in the possibility 
of computation of time and SM canonical averages of geometrical quantities entering the evolution of perturbations, 
used to determine the stability of trajectories avoiding the explicit integration of the variational equations, it emerges 
clearly why we need a recipe to cope with such singularities. This problem is also related to the Statistical Mechanical 



properties of gravitational N-body systems that we will not touch here, | Cipriani and Di Bari 1997b |. If we perform 
a static average of the Ricci curvature over M, we neither need nor can to avoid these 0{Af^) singularities, which 
are nonetheless Lebesgue-integrahle, and do not create any difficulty in the definition of the average values {kpi) m or 
any other related quantity; although they contribute heavily to the magnitude of fluctuations around these averages. 
Vice versa, when we look at time averages, we face up to a problem of probability: what is the chance that a true 
collision among particles do occur? If we have to do with the purely mathematical problem, we are without hope: the 
GDI's possess their well defined static averages, to which the corresponding time averages never approach! We feel to 
be faced again, from a novel perspective, to the classical problem of the non existence of an equilibrium state for the 
mathematical Newtonian gravity, dressed within this framework as the non existence of an ergodicity time. 

The way out to this problem is, as expected, a softening in the potential. This is an usual trick (either implicit or 
explicit) of N-body numerical simulators, adopted to cope with the so-called close encounters which occur sometimes 
during simulations. It is traditionally passed on that the occurrence of such kind of events is a consequence of the 
necessarily small number of particles used in a feasible simulation. This is true to a limited extent, from a numerical 
point of view, but it is a general occurrence, from a physical perspective, that during an orbital period some few stars 
suffer of such an event. We avoid to discuss here the details of the criteria adopted to define such a close encounter, 



[Cipriani and Di Bari 1997b |, and limit ourselves to say that when the distance between any pair of stars becomes 
exceedingly small with respect to the average interparticle separation, d — D/N^^^, the newtonian interaction 

should break down, as other effects, possibly even non conservative, take place. 

The easiest way to take into account for this breakdown consist to modify slightly (very very slightly) the interaction 
potential, allowing for a softening parameter, uj (we don't use the usual e to avoid confusion with the specific energy 
e) enter in the distance between particles used to compute the potential: 

rik =' [{xk - x^)^ + iVk - Utf + {zk - Zif + w"^] (61) 
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The introduction of the softening, helps to improve the conceptual understanding of the pecuharities of Gravity 
without modifying appreciably the qualitative overall dynamics, which is determined by the long range nature of the 
interaction, rather than by the short range divergencies, involving only a negligible fraction of masses (when 3> 1). 
The non-compact nature of the phase space, the bad statistical properties of gravitational interaction (that of being 
not a stable one) are indeed preserved by the introduction of the softening. This because it do not changes at all 
the long range behaviour of the interaction, neither introduces an hard core, as it could be erroneusly understood. 
The presence of the term lu^ in the potential has so theoretical consequences by far more relevant than the practical 
ones obtained in the more well behaved numerical integration of equations of motion. Indeed, this trick allows both 
to evaluate analitically and to numerically compute, just for a check on those estimates, everywhere non-divergent 
quantities entering the GDI's, to find their scaling law behaviour with relevant global parameter of the system (i.e., the 
number of degrees of freedom Af — 3N and the energy density e), by means of which to locate the possible relationships 
with the dynamical time-scale, t]j. 

As we will see, the analytical estimates suffice to single out what is it the dominant term in the frequency Q{t) appearing 
in equation (|3|), and also a precise determination of its relationship with the energy density e. For standard potentials, 
as the fP\J-(3 one, we have seen above that, provided Af ^ 1, neither the dynamical time to, nor the average values 
of kji and Q depends on A/", which weakly affects only the amplitude of the fluctuations. As expected, instead, the 
Gravity, ever corrected from too mathematical schematizations, shows a very different behaviour, disclosing a very 
intricate dependence of the GDI's on the number of particles N. This relationship can only be guessed through order 
of magnitude estimates, and we need to appeal to numerical simulations, whose results let us to state that: 

• The frequency Q{t) is overwhelmingly positive for any value of N, so the Chaos observed (or better, the expo- 
nential instability) is certainly due to parametric resonance. 

• The Ricci curvature fci^(s), and consequently k^it), is sign fluctuating for very small values of N, but shows a 
clear tendency to become more and more positive as the number of point masses increases. Already for some 
tens of bodies, it turns out to be mostly positive. 

• In any case, the evidence of the scaling law with N of all the terms appearing in kji{s), allows us to claim 
that, even for relatively small N, the average value of Ricci curvature is always positive, because the only term 
contributing with a constant negative sign is the faster to decrease with N. 

• The evidence of a virial transition with the implied damping of fluctuations. 

• If, moreover, the question of equivalence among static and dynamic averages is addressed, taking so into con- 
sideration the laplacian term, which is always the biggest one in the average, irrespective on the detailed way it 
is accounted for, the scaling law behaviour with N , prove that the Ricci curvature of the gravitational N-body 
system is almost everywhere (i.e., almost always) positive even for moderately large N. In the very limit iV ^ 1, 
kii{s) is always rigorously positive. 

• Well then, the geodesic instability too arise from the fluctuations of the positive Ricci curvature, which in 
turn implies the fluctuations (i.e., the anisotropy) of the sectional curvatures of the manifold. In such a case, 
neither rigorous results nor convincing arguments exist compelling the instability time-scale, of the order of the 
dynamical one, to give a guess on the SM relaxation time-scale. 

• The absence of any transition, analogous to the SST observed in the FPU model, in the degree of chaoticity of 
the dynamics, leaves open the question of the existence of different regimes of rapidity in the approach to the 
equilibrium. 

• The amplitude of fluctuations and the level of anisotropy in the iV-body system are nevertheless amplified with 
respect to those occurring in the FPU chain. It is not hard to guess for the sources of the strong anisotropy 
and of the big fluctuations in the configuration manifold. The results obtained for few dimensional DS's, 
[[Cipriani and Di Bari 1997c |, let us to support the claim according to which a self gravitating A^-body system 



find itself in a regime of fully developed chaos, for the vast majority of initial conditions, as long as TV ^ 1. 

The discussion reported above is the starting point for an investigation on the issue of the existence of one final 
or various partial SM equilibria for TV-body self gravitating systems. The analysis carried out on the relevance of 
the scaling law of the different terms appearing in the GDI's suggests a scenario with a full hierarchy of processes 
of approach to metastable equilibria, whose time scales are ordered in accordance to a criterion related to the 
scale (i.e., the fraction n/N of the stars involved) of the process itself, and will appear soon. 
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Comparison between analytical estimates and numerical results. 



In the table 1^ we report the analytical estimates, slightly improved with respect to those appearing in | Cipriani 1993 1, 



on the relevance and the dependence of the terms entering the GDI's, whose full account will appear in 

[ Cipriani and Di Bari 1997b| , along with some consistency checks obtained by numerical simulations. Minor cor- 



rections, which do not change the overall behaviours, have been obtained taking into account a more careful statistical 
evaluation of the amplitude of those terms we named fluctuating. As it is evident from the table, for the FPU-/3 chain 
there are no doubt on the relative weights of various terms. When the virialization phase is completed, and this happen 
at once, when the number of oscillators is appreciably greater than few ten^^ the average positive values of both Q{t) 
and kii{t) are determined mostly by the laplacian term, while the fluctuations around these values, again after the first 
phases of approach to virial equilibrium, are due exclusively to the laplacian for Q{t), whereas are influenced also by 
other terms for what concerns kR{t). Nevertheless, although characterized by greater fluctuations, the average values 
of ki{(t) obviously coincide with {Q(t)), as shown in fig.^. Moreover, as shown in fig.|^, after the virial equilibrium has 
been attained, as long as we have A/" 3> 1 and the system is far from integrability, the probability of the occurence of a 
negative value of the Ricci curvature becomes actually zero. Apart from this last remark, these results are substantially 
unaffected whether be the regime we are investigating, as all the terms scale with /3e in the same fashion, in such a 
way that their relative weights are unchanged. While, obviously, the amplitude of the fluctuations increases along 
with anharmonicity. 

So, for the FPUt3 chain there aren't many cautions to be made, and the numerical simulations simply confirm the ana- 
lytical estimations and support all the argumentations based on them, Cipriani and Pucacco 19944 Cipriani and Pucacco 1994b 



As a simple example we plot in figure ^, from [ Cipriani 1993| , the e dependences of (Q) and {ku), for three different 
values of /3, which exploits neatly the foreseen scaling as (/3e)^/^. 

As remarked above, the situation for the A^-body problem is more involved, and the analytical estimates alone let 
us know only some of the answers, although perhaps the most physically relevant ones, namely, where come from the 
most important contribution to Q(t) and the scaling behaviours of the Ricci related quantities with the energy density 
[e[, which lead to the conclusion that all them depend on the energetic contents of the system exactly as they have to 
do; for example, that the average value, measured in natural units, (Q) t^ do not depends on e. Stated otherwhise, 
this indicates that the scale invariance of gravity manifests itself within geometrical transcription in such a way that 
all the (possibly suitably rescaled to meet dimensional needs) averages of GDI's have the same values if measured in 
the unit of dynamical time scale t d . 

Moreover, having singled out the dominant contributions to Q{t), we are able to definitely support our claims that 
even in this case is the parametric resonance to drive to the onset of Chaos, and that the system dynamics should 
undergo a transition along the approach to the virial, if the initial conditions are chosen out of equilibrium. 
What it is left to determine is the scaling law behaviour of the GDI's with A^, which is practically absent in extensive 
potentials, but nevertheless likely to exist in the superextensive N-body problem, as indeed shown in table ^ The 
correct behaviour has been established supplementing the analytical investigations with numerical simulations whose 
details will be presented elsewhere, but whose main outcomes are: 



The frequency Q(t) is in almost always positive for the N-body problem, and the chances to find a negative value 
decrease quickly when TV increases. To show this, we report in figures |5[ ^ and |^ the time evolution of all the 
contributions coming from static and Ructuating terms. 

Despite the positive values of both Q{t) and fc_R(i) (whose average is always positive and for which also the 
occurrence of local negative values becomes negligibly small when A^ increases) imply that also the Ricci curvature 
per degree of freedom, fc/? (s), is positive; we found that the dynamics of self gravitating A^-body systems is always 



strongly chaotic, [Cipriani and Di Bari 1997b |, as far as chaoticity can be defined for non-compact phase spaces. 
This is explained again by the phenomenon of parametric resonance which is always in the regime of fully 
developed instability, as the fiuctuations in the GDI's are for this DS constantly above the threshold required 
for the mechanism to onset. 

All quantities entering the GDI's have the same specific energy e dependence as tj^^, i.e., both the laplacian and 
gradient related quantities, on one side, and the first (squared) and second time derivatives of the total potential 



We remind that the virial equations in their general form: {q°'dH / dq"") = {p'^dH/dp"'), where the averages are taken over time or 
phase space, when applied to mdf system, using the central limit theorem and summing over the degrees of freedom, let to state that 
the equality is true istantaneously with an approximation that becomes more and more reliable as N increases, being the statistical error 
proportional to N~^, where, depending on the interaction, ^ ^ b < 3, see also note ISl 
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energy, on the other, scale exactly as |ep, in such a way that the frequencies they define are in constant ratio 
with the natural orbital one. 



• The N dependence of these quantities, even when removed that coming from the dynamical time (we recall that 
we have to oc A^|e|~'^/^), manifests a very pecuhar behaviour, with all quantities changing with N, but in such 
a way that the hierarchy of their relative ratios either remains unchanged or becomes even more pronounced; 
figure js] shows the behaviour of the contributors to the Ricci related frequencies when N is varied. 

• With the support of specifically devoted numerical simulations, we explained this peculiar iV-dependence, which 
turns out to be very enlightening also on the consequences of the gravitational clustering, whose basic more 
elementary example is the formation of a binary system. The implications of this point on the evolutionary 
processes of stellar systems are currently under study. 

• The numerical simulations performed confirmed that the analytical estimates correctly describe the behaviour 
with N of the ratios between the terms containing the time-derivatives of the kinetic energy of the system. 

• Also for the static terms the a priori estimates correctly indicate their relative weights and even their ratios to 
the Ructuating ones. 



We investigated in detail also the behaviour of the static terms, and interpreted their behaviour with N on 
the basis of an analytical conjecture, confirmed by numerical computations, whose details will be discussed in 
[ [Cipriani and Di Bari 1997b |, where we will show also how all these outcomes do not depend on the detailed 



way we adopted to cope with the mathematical singularity related to the laplacian term in order to restore its 
physically relevant role. 

Figure ^ clearly shows the precise scaling law of all the quantities entering the GDI's on TV, along with the trivial 
dependence on |e|, or, is is the same, on ti). There are other minor comments which can be made to the results of 
numerical simulations, among which we mention the very small oscillations in the panels a) and c) of figure ^, along 
with their long overall quasi-periodicity (the time is in millions of harmonic units!). Moreover, panels b) and d) of 
the same figure point out how the relationship between Chaos and frequency of a negative sign of the Ricci curvature 
should be at least, bimodal. 



Conclusion. 



This paper belongs to a line of research addressed to the issue of a geometric description of Chaos in DS's. Here we 
focused our attention mainly of mdf hamiltonian systems, in order to investigate the sources of instability and to single 
out possible a priori criteria to detect the onset of chaotic behaviour, in analogy with what obtained in the case of 



few dimensional DS's, [ Cipriani and Di Bari 1997c |. 

We discussed how the criteria borrowed from abstract Ergodic theory could be used, if ever a realistic physical model 
would fulfil their rather restrictive hypotheses. To this goal we exhaustively investigated, starting at a very basic level, 
the relationships existing between the various curvatures that can be defined over a manifold M where the geodesic 
flow translating the dynamics in a geometrical picture takes place. We derived some of their logical consequences, and 
then indicated some cautious remarks in order to get reliable indications from ambitious applications of the GDA, 



[Gurzadyan and Savvidy 1986, Kandrup 1990 1. 

In particular, we here reported the formal explanation (see,e . g. , [Cipriani 1993|), within a more general context, of the 



analytical and numerical results we obtained, [Cipriani and Pucacco 1994a, Di Bari and Cipriani 1997c, Cipriani and Di Bari 1997c 



about the irrelevance of scalar curvature and related quantities to detect the presence of instability, obviously as long 
as A/" > 3. 

It is only recently, e.g., [ Pcttini 1993| , Cipriani 1993 [, that a careful framework has started to be set up. Thanks to 
this more systematic approach, the method has revealed, and we claim only in part, its fruitfulness, for both mdf and 
few dimensional DS's. 

In the present work we exploited its ability to single out the transition between two different regimes of Chaos in 
the fPi)^ model. Further, we also reinterpreted some sources of misunderstandings for what concerns the study of 
instability in the gravitational N-body problem. We shown its reliability performing analytical calculations, aided 
and supported by numerical computations, able to detect the scaling law behaviour of relevant quantities entering 
the equations governing the second order variational equations and in such a way identifying once more the peculiar 
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behaviour on Newtonian Gravity. We proved nevertheless that the mechanisms responsible of the onset of Chaos in 
standard mdf Hamitonian systems of interest in a classical Statistical Mechanics context, are also at work to determine 
the exponential instability (if not true Chaos), in this not well behaved DS. In this way, the somewhere persistent 
claims, |Gurzadyan and Savvidy 1986, Kandrup 1990|, about a naive relationship between instability and relaxation 
time scales should be revisited, and that their identification should be considered definitely false. Rather, it has to 
be reminded that for some DS's there are some analogies in their behaviours, though limited to the existence of a 
common transition between two different regimes, [Pettini and Cerruti-Sola 1991 , but that, for the N-body system 



no relationship has any more than heuristic support, [ Cipriani and Pucacco 1994a , [Cipriani and Pucacco 19941 1 



The GDA has made possible, the analytical evaluation of the maximal LCN for the FPU chain, through simple proce- 
dures developed in the framework of stochastic differential equations, [ Van Kampcn 19761 , using solely the time average 



of the Ricci curvature and of its dispersion along with an estimate of its autocorrelation time (see [Cipriani 1993 



Cipriani and Di Bari 1997b [ and [ Casetti and Pettini 1995 |). The analytically calculated maximal LCN's are in a def- 
initely satisfactory quantitative agreement with those numerically computed using the standard tangent dynamics 
equations, except at a very low (specific) energy, i.e., when Chaos is not yet fully developed, as in this case the Euctu- 
ations of GDI's show a slow convergence to their asymptotic values (whereas the average values of the curvatures are 
reached very quickly also in the harmonic regime) . Incidentally, this is a signature of the many faces of the approach to 
the equilibrium, which is a concept strongly dependent on the observable you are looking to relax to its ergodic value. 
Recalling a previous discussion, we furthermore observe that similar results are obtained averaging over phase space, 
instead of on time, and it should be intriguing to ask how the GDA allows also to gain some hints on the relaxation 
time scales for mdf systems whose final equilibrium is guaranteed to exist. 

The framework outlined here and previously developed in [Cipriani 1993, Di Bari 1996| serves as a starting point to 
try to answer also these basic questions. In the near future wc will exploit the results obtained both in the case of few 
dimensional DS's and for mdf lagrangian systems, to shed light, from a novel perspective, to the long standing issue 
of a dynamical justification of the statistical approach in (classical) mechanics. 
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Figure Captions 



Figure 1: Long time behaviour of the measured cumulative frequency of occurrence of negative values of the Ricci 
curvature F-{t) in the FPU-/3 chain, for a set of initial conditions far enough from integrability. All the solid lines 
refer to = 450 coupled oscillators with fie ranging from 5 • 10~^ to 10'^. The dashed lines represent the behaviour of 
a smaller system {N = 150) not too far from integrability, j3e = 10"'^, lO^'^, from bottom to top, showing the greater 
noise affecting simulations with too small N . 

Figure 2: The same of the figure |l|, but in the opposite situation: very small N and/or conditions very near to 
integrability. The values of {N,/3e) are: a) (50,5 • 10~^); b) (50,0.03) (upper curve) and (50,0.01) (lower curve); c) 
(50,5 • 10""*); d), from top to bottom, (150,5 • IG"'^)^ (^50, lO"'^), (450, 5 • 10"^) and (150,2 • lO"'^). The horizontal 
scale now is always linear and time is there expressed in million of units. 

Figure 3: Energy dependence of the largest LCN of the fP\J-(] chain, computed in analytical way. The dashed lines 
trace the low and high energy slopes: (/?e)^ and (/3e)^/'', respectively. In the small panel we plotted the asymptotic 
values of the cumulative frequency F_ {t) for the same conditions used to compute the LCN's. 

Figure 4: Time averages of the rescaled Ricci curvature k]i(t) and of the Hill frequency Q{t) for the same set of initial 
conditions use to compute the largest LCN. The continuous line is the plot of the analytic estimate reported in table 
The SST is neatly visible. 

Figure 5: Time behaviour of the quantitities entering the GDI's, Q{t) or kfi{t), for the gravitational N-body system 
with N = 10 (!) and |e| = 1. The small panels show, from left to right, the contributions to Q{t) (shown alone in 
the larger window) from the fluctuating and static terms, respectively. In the small panels here and in the following 
figures, the laplacian term is indicated with a solid line, the squared gradient as dot-dashed, the (positive values) of 
the second derivative of the potential energy as dashed, and the dotted curve display its squared first time derivative. 
In this and the following figures all the terms are plotted with the weights they enter in the expression of Q{t). 

Figure 6: The same as figure ||, but for N ~ 100. Note the logaritmic vertical scale of the inserted plot. The quantities 
are represented as in the previous figure. 

Figure 7: The same as the previous two figures but for N — 200 and with all the quantities in a plot. The small 
solid triangles trace the behaviour of Q{t) in order to try to disentangle it from that of ALl /A/", with which practically 
coincides. The other lines confirm the relative hierarchy of other terms. 

Figure 8: In the large area we plot the A^ dependence of the time averages of the various terms. These different scaling 
with A^ are explictly shown in the inserted plot, which also points out the sharp |ep-dependence of all these quantities. 



28 



Tables 



Table 1: Analytical estimates of the {N, e)-depenclence of quantities entering the GDI's and of the dynamical time, for 
the FPU-/? chain and the Gravitational N-body system, along with the indication about the scaling law behaviours 
observed for the instability exponents of both dynamical systems, evaluated either numerically or using the semi- 
analytical method discussed in the text. In this table N represents the number of oscillators or the number of masses. 
We remark that the amplitudes here indicated should be taken as order of magnitude estimates, up to numerical 
factors of order unity. Moreover, for what concerns the fluctuating terms, these evaluations represent almost the 
largest absolute value, as they average can be much smaller (even vanishing!). The quantities indicated with the 
symbols Cf„ depend on the clustering level of the N-body system and then slowly evolve, is a tracer of the way we 
adopted to cope with the newtonian singularity and the exponents a„ reflect the efiiciency of the virial damping on 
the amplitudes of evolving terms. Their explicit expression, as well that of the analytical formula, A, to compute the 
LCN, is given in jCipriani 1993, Cipriani and Di Bari 1997b|, see also [Casetti and Pettini 1995{. 
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